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During  project  PARKA,  data  -were  received,  digitized  and  processed 
in  real  time  by  the  UNIVAC  1230  computer  aboard  the  prime  receiving  vessel 
SANDS.  The  final  form  of  the  processed  data  was  displayed  as  listings 
and  graphs  with  the  values  being  expressed  in  /lb.  To  convert  the  data 
to  db,  two  generalized  programs,  ^TENLOGTEN*'*^ and  *C0NL0GIT*  were  written 
to  calculate  P = 10  log^o  (X) , where  X is  a fixed  point  number  within  a 
specified  set  of  limits.  The  major  differences  between  these  two  routines 
are  the  methods  used  to  obtain  the  natural  logarithm  and  the  scaling 
necessary  for  the  input  values. 

\ 

ADMINISTRATIVE  INFORMATION 

This  memorandum  was  prepared  under  NUSC  Project  Title:  Long-Range 
Acoustic  Transmission  Experiments  for  Surveillance  Systems  Development; 

R.  Hasse  and  R.  Martin,  NUSC/NL  Principal  Investigators.  The  sponsoring 
activity  was  0NR,  Code  102-0S;  Dr.  J.  B.  Hersey,  Program  Manager. 

"TENL0GTEN" 

Routine  "TENLOOTEN"  calculates  10  log^o  X by  evaluating  a 5th  degree 
polynomial.  Expressing  10  logiO  X as 


* This  routine  was  written  for  NUSL  by  TRAC0R  under  contract  for  project 
PARKA.  
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= (10  log10e)as  = .13966237. 

Routine  "TENLOGTEN"  calculates  N and  evaluates  the  above  equation. 
A flow  chart  of  "TENLOGTEN"  is  shown  in  Figure  1. 


Program  Usage 

To  call  on  routine  "TENLOGTEN",  enter  the  value  X in  the  A register, 
scaled  22  bits,  and  RJP*TENLOGTEN.  The  answer  appears  in  the  Q register 
scaled  3 bits. 

If  the  inputted  number  X is  less  than  or  equal  to  zero,  the  routine 
returns  a value  of  4000000000g. 


"CONLOGIT" 


"CONLOGIT"  also  calculates  P = 10  log^o(X),  by  obtaining  the  nat- 
ural logarithm  using  a subroutine  called  "NATLOG"*  and  converting  to 
P by  equation  (3).  To  use  "NATLOG",  the  number  X must  be  expressed 
in  the  form  X = Nowhere  l<N£377777g.  "CONLOGIT",  assuming  the 
scale  factor,  S is  15,  checks  N to  see  that  it  is  within  limits.  If  it 
is  not,  "CONLOGIT"  shifts  N down  so  that  the  most  significant  bit  is  in 
the  proper  position  and  then  adjusts  S accordingly. 


"NATLOG"  uses  a method  attributed  to  Dr.  H.  Maehly**  to  calculate 
the  natural  logarithm.  Rewriting  X as 


X = 2m  F, 


(5) 


where  m is  an  integer  and  li  Fi2, 
ln(X)  = (m+-2’)  ln2  + Po  + 
where  the  vertical  bars  imply  'continued  fractions,'  that  is 


(6) 


t 


- P, 


r«T. 


r * 


(7) 


* NATLOG  was  supplied  to  NUSC  by  UNIVAX  as  part  of  the  software  package 
delivered  with  the  1230  computer. 

**  H.  Maehly,  Monthly  Report,  Institute  for  Advanced  Study,  Princeton, 
Oct  1956. 
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and  Pq  = 3*681656603, 

P = 34.U1069291, 

P2  = 8. 12650383^, 

P3  = 0.2665195666, 

Tj.  = 10.37967214, 

T0  = 2.051212813,  and 
T3  = 0.4249952497. 

Figure  2 is  a flow  diagram  of  "CONLOGIT"  and  Figure  3 is  the  flow 
chart  of  "NATLOG". 

Program  Usage 

To  use  "CONLOGIT",  the  number  X is  put  in  the  Q,  register, 
scaled  to  15  bits.  If  X is  negative,  "CONLOGIT"  sets  it  positive  and 
proceeds,  or  if  X is  zero,  an  answer  of  -50.0  is  given.  If  an  error 
exit  occurs  in  "NATLOG"  caused  by  the  scale  factor  S being  outside  of 
the  limits,  an  answer  of  -100.0  is  given.  The  answer,  10  log10(x)  will 
be  in  the  Q register  scaled  3 Mts. 

To  call  on  "CONLOGIT",  use  the  call  "CONLOGIT"  or  use  "RJP*C0NL0' 


CLAIR  J.  BECKS6 
Mathematician 


DAVID  M.  POTTER 
Mathematician 
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APPENDIX  A 

LISTING  OF  "TENLOGTEN" 


TENLOGTEN 

ENTRY 

INPUT  IN  A SCALED  22D 

ENT#  B7  » LOG  CO  UNT 

N GOES  TO  W( LOGCOUNT)  ON  NORM  INSTRUCTION 

ENT#Q«0#AP0S 

PREPARE  Q FOR  NORMALIZATION 

JP# LOG ERR 

NEGATIVE  INPUT 

JP*  LOG  ERR*  AZERO 

INPUT  WAS  0 

NORM# AQ#W(B7) 

1 . GE  . X . LT  . 2 

SUB# A# 2000000000 

STR# A# W (LOGARG) 
ENT#B7#3 

MINUS  1 

0 

ENT#Q#W(L0GC0NT4) 

B4 

• 

LOGLOOP 

MUL#W (LOGARG) 

EVALUATE  POLYNOMIAL 

A 

ADD#  A#W(L0GC0NTB7) 

ADD  COEFFICIENTS 

STR#A#Q 

SCALING  IS  27,25,23,21  AS 

BJP#B7#L0GL00P 

B7  IS  3, 2, 1,0 

MUL#W(LOGARG) 

S TRV A *W( LOGARG) 

POLYNOMIAL 

ENT#  Q#  W ( LOG  CO  UNT ) 

N 

SUB#Q#6 

N-6 

MUL#W(L0GC0NT5) 

(N-6) (-3.0102) 

ADD# Q#W (LOGARG )*QPOS 

PLUS  POLYNOMIAL 

SUB#Q#100000#SK1P 

ROUND  OFF 

ADD# Q# 100000 
RSH#Q#16D*SKIP 

SCALE  TO  3 PLACES 

LOGERR 

ENT# Q# 4000000000 

ERROR  VALUE 

EXIT 

OUTPUT  IN  Q SCALED  3 

LOGCON 

DATA# 4 . 3407541D,21D 

DATA# -2  . 1363335D,23D 
DATA«1  . 2571730D,25D 
DATA#-. 59091 302D,27D 
DATA#. 139662 37D,29D 
DATA# -3  . 0102D,19D 
LOGARG  0 

LOGCOUNT  0 
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APPENDIX  B 

LISTING  OF  "CONLOGIT" 


CL2 

PROCEDURE* CONLOGIT 
3TR#Q*  A*  APOS 

CP*  Q 

JP*CL4*AZER0 

STR*B7*L(CL7) 

ENT*B7#1 

LSH*Q*1 

CL3 

JP# CL3*QNEG 
BSK*B7#13D 

JP*CL2 

JP#CL5 

ENT*B7*B7-13D 

CL  5 

STR*B7*CPL(l67) 

ENT*Q*X(B7-15D) 

ai 

CP*Q 

RSH* A*B7 

RJP* NATLOG 

CL  6 

JP*  CL6 

RSH*AQ#30D 

MUL* 33626754 
RSH*AQ#3bD 

MUL#12 

ADD*Q* 400* APOS 
SUB*Q*  4000 
RSH*Q#yD*SKIP 
ENT»Q*X76337 

CL7 

ENT*B7*0 

CL4 

RETURN 

ENT»Q*X77157 

JP*L( CONLOGIT) 
END-PROC#  CONLOG IT 

IS  NUMBER  NEGATIVE 
YES,  COMPLEMENT 
IF  ZERO,  JUMP  TO  CL 4 
SAVE  B7 

INITIALIZE  COUNTER 
CHECK  FOR  POSSIBLE  SCALE  DOWN 
SCALE  DOWN  REQUIRED,  JUMP  TO  CL3 
ALL  0\  RFLOW  VALUES  CHECKED? 

NO 

NUMBER  WITHIN  RANGE 
SET  SCALE  FACTOR 
MAKE  POSITIVE 
SET  NEW  SCALE  FACTOR 
MAKE  POSITIVE 
RESET  INPUT  VALUE  IN  A 
COMPUTE  NATURAL  LOG 
ERROR  EXiT  FROM  NATLOG 
NAT  LOG  IN  Q SC=24D  BITS 
CONVERT  TO  COM  LOG  SC=18D 
COM  LOG  IN  Q SC=12D  BITS 
CONVERT  TO  10  log  (db) 

ROUND  OFF  POSITIVE  NUMBER 
ROUND  OFF  NEGATIVE  NUMBER 
SC=3  BITS 

ERROR,  SET  ANS  TO  -100. Odb 

RESTORE  B7 

EXIT 

-50.0  IF  NlMBER-0 
EXIT 
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NATLOG 


APPENDIX  C 


LISTING  OF  "NATLOG" 


JP*0 

ENTRANCE 

STR*  Q* CP W (KITTY ) 

COMPLEMENT,  SAVE  SCALE  FACTOR 

CL*Q 

CLEAR  REGISTER 

ENT*  A*  A*  APOS 

IS  NUMBER  POSITIVE 

CP*  A 

RPT*29D 

NO,  COMPLEMENT 

LSH*A*1«  ANEG 

AND  COUNT  SHIFT 

JP*NAT2 

NUMBER  IS  ZERO 

STR*B7*L (NATl) 

LSH*  A*  W (K  ITT  Y+2 ) 
ENT*Q*X(0) 

SAVE  COUNT 

ADD*Q*W(KITTY) 

ORIGINAL  INPUT 

LSH*Q*3 

M SC  3 IN  Q 

ADD*Q*4 

CL*  A 

M+£  SC=3 

MUL*W(POOL) 

(M+£)LN2  SC=3+29D 

RSH*AQ#9D 

ENT*  A*  A*  AZ  ERO 

RIGHT  TO  SC-23D 

CP*A*ANOT 

IS  A 0 or  all  7's 

RPL*Y+1*L (NATLOG) *SKIP 

YES 

JP*L (NATLOG) 

STR* Q*W(KITTY+l) 

ERROR  EXIT 

ENT*Q*W(KITTY+2) 

F 

ENT*Y+Q*W(P00L3) 

F+T3  SC=23D 

ENT* Y+Q*W(P00L3) 

STR* A*W(KITTY+3) 

F+T3  SC=23D 

ENT* Y+Q*W(P00L2) 

STR* A*W(KITTY+4) 

F+T2  SC=23D 

ENT*  Y+Q*  W (P00L1 ) 

STR* A*W(KITTY+S) 

F+Tl  SC=23D 

ENT*  Q*  W (POOL  3+1 ) 

-P3  SC=23D 

EOT*  A*X77777 

SET  SIGN 

LSI!*  AQ*23D 

DIV*W(KITTY+3) 

ADJUST  SC  TO  4fD 

RPL*7+Q*W(KITTY+4) 

SAVE  FOR  NEXT  CONTINUED  FRACTION 

ENT*C,W(P00L2+l) 

-P2  SC-2 3D 

ENT* A*X77777 

SET  SIGN 

LSH*AQ*23D 

DIV*W(KITTY+4) 

ADJUST  SC  TO  46D 

RTL;-  Y+Q»W(KITTY+3) 

SAVE  FOR  NEXT  CONTINUED  FRACTION 

ENT»y»W(POOLl-*l) 

-Pi  SC=23D 

ENT»A»X77777 

SET  SIGN 

LSH*AQ»?3D 

DIV*W(KITTY+'5) 

ADD*Q*W(KITTY+l) 

ADJUST  SC  TO  A6D 

ENT*  Y Hi*  W (POOI.+ : ) 

ANSWER 

LSH»AM 

SCAI.E  ANSWER  TO  2.J)  BITS 

APPENDIX  0 (cont) 
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JP  *L ( NATLOG ) 

NORMAL  EXIT 

POOL 

2613441377 

LN  (2 ) SC=29D 

0166640206 

PO  SC=23D 

POOL  3 

0015463077 

T3  SC=23D 

7767361257 

-P3  SC=23i> 

P00L2 

0101507044 

T2  SC=23D 

7373747270 

-P?  SC-23D 

P00L1 

0514114431 

T1  3C-23D 

5662667151 

-PI  S 0=2 3D 

KITTY 

RESERVE  6 
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